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Abstract: We give a short overview of recent developments in understanding of the 
deeply virtual Compton scattering on the proton target. 
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Deeply virtual Compton scattering (DVCS), recently measured by HERMES Q, HI Q and C ^ 

CLAS Q collaborations at DESY and Jefferson Lab, is the cleanest hadronic reaction that ^ ^ 

gives access to the generalized parton distributions (GPDs) 0|. GPDs are the probability 
amplitudes to knock out a parton from, say, the nucleon and to put it back with a different 
longitudinal momentum fraction. The probed partonic state is characterized, of course, 
by quantum numbers, like flavor, spin, etc. GPDs unify a number of known concepts of 
hadronic physics: they are related to parton densities, wave functions, and parton form 
factors, coupled also to higher spin probes. Recall that conventional form factors are 
encoded in matrix elements of local currents without derivatives, while for the case at 
hand, operators contain any number n of derivatives. For n = 1 the form factors carry 
information on the orbital angular momentum carried by constituents in the proton, see the 
second paper of Ref. [Q]. Apart from DVCS and other related two-photon processes, GPDs 
also contribute to the hard leptoproduction of mesons, however, their theoretical description 
involves a new unknown nonperturbative input, a meson distributions amplitude, which 
complicates the disentanglement of GPDs from measurements. 



2. Differential cross section of eA^ e'N'j 

Presently, we concentrate on the structure of the cross section for electroproduction of the 

real photon off the proton e{k)N[Pi) — > e{k')N{P2)j{q2), whose amplitude is the sum of 
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Bethe-Heitler (BH) and the wanted DVCS amplitude. The four-fold cross sectfon 
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depends on the Bjorken variable xb = Q^/(2Pi • qi), the squared t-channel momentum 
transfer = {P2 — Pi)^, the lepton energy loss y = Pi ■ qi/ Pi ■ k, and the azimuthal angle 
(j), see Fig. ||. The resolution scale = —qf is given by the virtuality of the incoming 
photon. In our frame, the angular dependence of separate components in the cross section 
(|2.1|) is given by a finite sum of Fourier harmonics: 
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Figure 1: The kinematics for leptoproduction 
in the target rest frame, where the azimuthal 
of the proton momentum with respect to the 
lepton scattering plane is cf). 



Here the + (— ) sign stands for electron (po- 
sitron) beam. In the squared BH and the 
interference term an additional (f) dependence 
arises from the scaled lepton BH propagators 
Vi{4>) = Ai + i3jCos((/>), which can be rather 
strong for large y. The coefficients Cn and Sn 
depend on the kinematical variables y, xb, 
A^, Q^, and lepton and hadron polarizations. 

The Fourier coefficients (FCs) cf and sf 
[cj, cf, and sf] as well as Cq^'^^ [c^^"-^^, and 
^pvcSj appeal- at the twist-two [-three] level, 
while even higher harmonics are suppressed 



by powers of Og or For an unpolarized target we have in the twist-two sector 
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(2.4) 



where K ~ y^(1 — xb)(1 — y)A_L / Q and A is the lepton helicity. The Cs are functions of 
the so-called Compton form factors (CFFs) J- = {Tl, £, T~L, £} which parametrize the DVCS 
tensor and they have been worked out in Refs. (H at twist-two level, e.g.. 
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where Fi and F2 are the nucleon Dirac and Pauh form factors, respectively. The twist-three 
analysis |^] (see also 0) has just been completed and will be presented in detail elsewhere. 
The CFFs are expressed in terms of GPDs F = {H,E,H, E}. At leading order (LO) of 
perturbation theory the relations read 

(2.5) 

where Qi is the fractional quark charge and ^ = xb/C^ — xb)- Analogous formulae are 
valid for the spin- flip CFFs S and £. Unfortunately, from Eq. (|2.5|) one can not practically 



deconvolute GPDs. However, as Eq. (2.4) demonstrates, those observables which are sen- 
sitive to allow to measure directly GPDs (at LO) for x = Generally, one has to rely 
on models with a set of free parameters, which have to be adjusted to experimental data. 

3. Observables and extraction of model parameters 

For the low-Q^ kinematical settings of the present experiments ||l|, ^] it is desirable to 
separate the twist-two and -tree sector, since a priori the latter can contaminate the former. 
This procedure is practically possible only with facilities where positive and negative lepton 
beams are available by means of charge asymmetries and Fourier analysis. 

, , „ „ The charge-odd combinations of cross sec- 

tion (|2.l| ) gives the interference term whose co- 
efficients can be extracted by the weighting pro- 
cedure. 
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Figure 2: Alu as a function of (j>'^ = 
TT — (/) at HERMES versus model predictions: 
AvKW (solid) and Jiww+qGq (dash-dotted). 
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(3.1) 

with the measure dw oc Vi{(p)V2{(p)d(j). The 
FCs of the squared DVCS term can be obtained 
from the charge-even part after subtraction of 
the squared BH term, i.e. d^^^^o" = {da~^ + 
d~a)/2 — d^^a, and weighting with respect to 
the measure d(j). Asymmetries are normalized 
which does not contain twist-three corrections. Unfortu- 



with respect to Jq 

nately, for single-charge lepton beam machines this clear separation of twist-two and -three 
effects cannot be achieved - also not with spin or azimuthal asymmetries. For instance, 
for the beam-spin asymmetry measured with HERMES and CLAS 

^Lu(0) = (da^ - da^^ / (^da^ + da^^ , (3.2) 

the leading twist term in the numerator will be affected by S2 ^^'^ s^^^^ while the 



denominator will be the sum of all three contributions in Eq. (2^). Fortunately, if the 
condition (1 — y)A^/y^Q^ <C 1 is fulfilled, s^^^^ is kinematically suppressed and, moreover, 
the squared BH term Cq ^ dominates in the denominator. Thus, the lepton BH propagators 
approximately cancel in Eq. ( ^I^ ) and ^lu(0) ~ i(a^BSi'/yco^) sin(0). 
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Since the number of experimental observables is very limited we are confined to use 
GPD models in order to confront theoretical expectations with data. Let us briefly outline, 
however, a way how model parameters can be fit to experiments in particular kinematical 
situations. An ansatz is based on a simplistic factorization of and (x, ^) dependence. 
Sum rules and reduction formulae and ad hoc assumptions result into a model for, e.g.. Hi 
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Qi are the parton densities given at the input scale Qo and N{hi) = -^^f^^^- 
For valence quarks the form factors are fixed by the sum rules, while for sea quarks we 
parametrize them as Ff^^{A.'^) = (l — i?sca^^/3) ^ with a free slope parameter i?sca- The 
free parameters bi in Eq. ( |3.3D control the skewedness effects. The last term in ( |3.3| ) stands 
for the mesonic-like contributions, see |^. Parametrizations for the spin-flip GPDs are 
obtained by the replacements: Fl — > F2 and Di —Di. Hi is given in terms of polarized 
parton densities, while Ei is dominated by the pion pole, see 

Assuming the MRS A' and GS A densities 
at the input scale Qq = 4 GeV^ , we analyzed the 
existing DVCS data with these models at LO 
in the twist-three approximation with neglected 
Q^-evolution. The unpolarized cross section at 
small xb as measured by HI Q is dominated by 

T^sea, i.e., CDVCS ^ |^|2 foj. _^2 ^ 4^j2_ ^e 

found that the skewedness effect is small and the 
fall-off with is strong. The choice 6sea — > 00 
and Bsca. = 9 GeV^ gives a good agreement with 
the data. Using 6vai and sea-quark magnetic 
moment Ksea as free parameters we distinguish 
Vai = 1 and Ksca = 0, and B with 6vai 00 and Kgea. = —3. 
In Fig. and |3| we give our predictions for the model A with twist-three GPDs taken 
in the Wandzura-Wilczek approximation (solid line) @ and the model B with quark- 
gluon-quark correlations accounted for in the twist-three GPDs (dash-dotted line) versus 
experimental data. For HERMES and CLAS plots, the dotted line denotes the sm(j) with 
the amplitudes 0.23 and 0.202, respectively, fit to the data by the collaborations. At 
HERMES the integrated beam-spin asymmetry Alu = -0.23 ± 0.04(stat) ± 0.03(syst) has 
been measured at the average values (Q2) = 2.6 GeV^, {xb) = 0.11 and (-A^) = 0.27 GeV^ 
Q . Integrating our predictions Avi/vf (solid) and Bww+qGq (dash-dotted) from Fig. |2| we 
have ^LU = —0.27 and Alu = —0.16, respectively. Both of them are compatible with the 
data. The beam-spin asymmetry as measured by CLAS collaboration was integrated over 
the region: 1 GeV^ < < 1.75 GeV^, 0.13 < xb < 0.35, and 0.1 GeV^ < -A^ < 0.3 GeV^ 
with W > 2 GeV. Their integrated result is fit by Ai^u{(f>) = asm{(j)) + (3sm{2(f)) with 
a = 0.202 ± 0.028(stat) ± 0.013(syst) and p = -0.024 ± 0.021(stat) ± 0.009(syst) [|. The 
model Av(/vF fails to describe the data resulting into a = 0.28 and /3 = 0.014. However for 
the model B^yy^gGq we find a = 0.24 and /? = —0.03 in fairly good agreement with the 



Figure 3: Alu as a function oi (p'^ = it — 
(j) at CLAS versus model predictions. Same 
assignments as in Fig. ||. 

two models: A with h 
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measurement. Note that contrary to a, the second harmonic parametrized by /3 is very 
sensitive to quark-gluon-quark correlations. 



4. Beyond leading order and power 

So far we have shown that the DVCS results measured in three different experiments are 
consistent for a set of GPD models with DVCS amplitudes evaluated at LO within the 
twist-three approximation. Once we will have a rough idea on the magnitude of model 
parameters a refined analysis has to be conducted taking into account higher order effects 
in QCD coupling constant. In the twist-two sector, they are completely worked out at NLO 
Q. According to our recent analysis, their size is moderate for each quark species, however, 
large radiative corrections can be generated by certain models of gluon GPDs. The latter 
are induced by strong evolution effects and, thus, these models can be experimentally tested 
in the small-ZB region. The gluonic contributions can be set to zero by an appropriate 
choice of the factorization scale. In this scheme we have moderate 20% corrections to the 
CFFs, however, the DVCS cross section then gets corrected by 40%. 

As we have observed the CLAS result presumably indicates a non-negligible contri- 
bution from multiparton correlations in the nucleon. Therefore, a further study of higher 
twist effects is necessary. Since is of the order of one GeV^, one expects also a sizable 
contamination by target mass M^/Q^ effects. A first step in this direction has been made 
in [^, however, a complete treatment requires a consideration of mass effects stemming 
form multiparticle operators. 
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